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The problem of diffusive bond- dissociation in a double well potential under application of 
an external force is scrutinized. We compute the probability distribution of rupture forces 
and present a detailed discussion of the influence of finite rebinding probabilities on the 
dynamic force spectrum. In particular, we focus on barrier crossing upon extension, i.e. 
under linearly increased load, and upon relaxation starting from completely separated bonds. 
For large loading rates the rupture force and the rejoining force depend on the loading rate 
in the expected manner determined by the shape of the potential. For small loading rates the 
mean forces obtained from pull and relax modes approach each other as the system reaches 
equilibrium. We investigate the dependence of the rupture force distributions and mean 
rupture forces on external parameters like cantilever stiffness and influence of a soft linker. 
We find that depending on the implementation of a soft linker the equilibrium rupture 
force is either unaffected by the presence of the linker or changes in a predictable way with 
the linker-compliance. Additionally, we show that it is possible to extract the equilibrium 
constant of the on- and off-rates from the determination of the equilibrium rupture forces. 



PACS numbers: 82.37.Np, 82.37.Rs, 87.10.Mn, 87.15.Fh 



I. Introduction 

The determination of rupture forces, e.g. via atomic force microscopy (AFM) or optical 
tweezers by now is a standard method of investigating the strength of molecular bonds, for 
reviews see[TJ [21 El II]- In ah the methods at hand one studies the breaking of a bond under 
the influence of an external force, mostly by applying a linear force ramp. The observed 
rupture forces are then monitored on a single molecule level. Therefore, one observes a 
distribution of rupture forces subject to further analysis. In most examples considered 
experimentally and theoretically bond rupture is treated as an irreversible event. This 
means that after a rupture event the molecular configuration or structure has changed 
irreversibly into a very different arrangement void of memory of the prior structure. The 
analysis of either the mean rupture forces or the distribution of rupture forces yields valu- 
able information about the details of the bond under study. Most theoretical models that 
have been developed to treat this dissociation process consider it as being a problem of 
diffusive barrier crossing[5]. Many studies have dealt with the detailed dependence of the 
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transition rates, i.e. the off-rates, on the applied force[6j [7J El E], on the elasticity of soft 
linkers [T0| ITT| [T2] and on intrinsic bond heterogeneities [13]. 

Although most biomolecules unfold irreversibly under the influence of mechanical force 
there are some important exceptions. Among these is the unfolding of RNA hairpins [TU [15] 
and of ubiquitin[16j and the unzipping of coiled coil peptides]!?]. The reversible bind- 
ing/unbinding of adhesion molecules has been studied only in recent years [TBI [19]. Seifert, 
using Kramers theory, has shown that rebinding has a strong impact on the rupture force 
versus loading rate relation and that equilibrium is asymptotically reached for small load- 
ing rates [TS]. Only in the limit of large loading rates rebinding effects become negligible 
and the behavior expected for purely irreversible rupture is recovered. Li and Leckband[T9"] 
expanded the analysis by a detailed discussion of the dependence of the calculated rupture 
forces on kinetic parameters for one specific model potential of a double well type. Both 
approaches treat the case of a large number of equivalent adhesion molecules. Therefore, 
the force-versus-extension (FE) curves or the so-called dynamic strength corresponds to 
'bulk' properties and the difference to properties observed in single molecule experiments 
has been mentioned already by Seifert [18]. A finite number of reversibly breaking adhesion 
bonds has been treated in a one-step master equation approach by Erdmann et al.[20]. 
These authors presented a detailed discussion of the corresponding rupture force distribu- 
tion. The analysis has some similarity with the one for the dissociation of sequences of 
identical proteins [8] with the difference of finite rebinding rates. 

In the present paper we consider the reversible rupture of single molecules that exists 
in two states, e.g. folded or unfolded. We treat the problem of diffusive barrier cross- 
ing in a double well potential similar to earlier studies. We compare the results of our 
semi-analytical calculations based on the adiabatic approximation to Brownian dynamics 
simulations. In contrast to existing studies we primarily focus on the distribution of rup- 
ture forces as obtained e.g. in AFM experiments and explicitly consider the behavior of 
these distributions in the situation where the external force is linearly reduced. We discuss 
the dependence of the resulting distributions and rupture /rejoining forces on the model 
parameters and also consider the influence of flexible linkers and bond heterogeneities. The 
remainder of this paper is organized as follows. In the next section we present the theo- 
retical calculations and section III. contains the discussion of the results. The paper closes 
with some concluding remarks putting emphasis in how data analysis should be carried 
out to avoid possible misinterpretations. 

II. Theory 

We consider reversible bond breaking as a process of diffusive barrier crossing. In the 
'standard' model one considers Brownian motion in a potential of the general form: 

V(q, t) = V (q) + V m n(q, t) with V m n(q, t) = ~k c (q - vtf (1) 
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Here, q denotes the reaction coordinate, k c is the force constant of the cantilever and v is 
the drift velocity. In the simplest case the force acting on the system is given by / = k c -v-t, 
typically in the pN- range in AFM experiments. The potential Vo(q) is chosen to be of a 
double well form with one minimum located at q% a transition state at and another 
minimum at q B (q B > q? > q A )i cf. Fig. la. In order to study the Brownian dynamics in 
such a potential, one has to solve the overdamped Langevin equation associated with the 
potential (JT|, given by: 

q(t) = - 1 V\q ) t)+0t) (2) 
where 7 is the damping constant and £(£) is delta-correlated Gaussian white noise, 

(£(*)) = an d (£(*)£(*')) = 2 7 T5(t-0 

This Gaussian property of the noise is the reason for the fact that the stochastic process 
q(t) is a Markov process. Equivalently, the process can be described by the Fokker-Planck 
equation (Smoluchowski equation) [21| [22] : 



p(q,t)=D 



^_ e ~(3V(g,t)^_ e l3V(g,t) 

dq dq 



p(q,t) (3) 



Here, (3 = 1/T and the diffusion constant is related to the damping constant via D = 7T 
(Boltzmanns constant is set to unity throughout). The solution of either equation has to be 
performed numerically, as no analytical solutions are available for the nonlinear potentials 
considered here. In this paper we will present results of Brownian dynamics simulations, 
i.e. the numerical solution of the Langevin equation, eq. (J2| [23J . 

From the Fokker-Planck equation one can derive expressions for the transition rates 
from A to T and from T to B as mean first passage times (MFPTs) [ 



tatU) = ^ ["dqe^ (* dq'e-^ ; r BT (f) = ± [** dqe*™ Hdq'e-^ (4) 

where all quantities qA, qr, Qb and V(q) are functions of the force /. If the potential barrier 
is not too low, it is a good approximation for tat{J) to replace the upper limit of the inner 
integral by the position of the transition state, qr(f)- For r BT (f), one replaces the lower 
limit by qr(f)- This way one obtains the product approximation: 

r P-duct (/) ^ D -i r qT dqe mi) [ qT dq'e-^ 

J qA J —00 

r P^ duct (y) ~ D~ x f qB dqe^ Hdq'e-^ (5) 

Jot J 



The transition rates for passing from potential well A to well B can be calculated from 
the MFPTs given above using a three state model A «-> T «-> B and assuming that the 
concentration of T is stationary. The effective transition rates are given by [5]: 

Ks(f) = TAT(f) + r BT (f)Z A (f)/Z B (f) 

KM = TBT(f) + r AT (f)Z B (f)/Z A (f) (6) 
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One particular example of a double well potential that allows the analytical calculation of 
the MFPTs according to eq.([5]) is given by a harmonic cusp-like potential, cf. Fig. la: 



v (n \ _ / V A + \k A {q ~ q\Y for q < q° T 
V ° [q) ~\V B + \k B {q - q B f for q > q° T 



(7) 



A similar potential has been considered by Li and Leckband in their discussion of the 
effect of rebinding on the average dynamical strengths of molecularly bonded surfaces [P9"]. 
In eq.([7]), the force constants are denoted by k^ and k B for the left (A) and the right 
(B) well, respectively. The q\, X = A, T, B are the positions of the extrema in the 
absence of force. Finally, the barrier is determined by the intersection of the two parabolas, 
Vt = l/2kA(qx) 2 • Note that in Fig. la the overall potential according to eq.((T]) is plotted. 
The explicit expressions for the MFPTs for the cusp-like potential given in eq.([7]) are 
provided in Appendix A. There also the Kramers rates for a general double well potential 
are given along with the limiting case, the rates in the Bell model[6j. 

In Fig. lb we show the rates k g(f) and k on (f) resulting from the exact expression, 
eq.Q, from the product approximation ^ and the Kramers approximation, eq.(A.4). It 



is evident that the product approximation yields excellent results, whereas the Kramers 
approximation fails when the force at which well A vanishes is approached. 

Denoting the time- dependent population of well A by n(t), then the master equation 
reads after a transformation to the force: 

^n(f) = - [n ofi (f) + K on (f)} n(f) + K OQ (f) (8) 

Here, we defined n n/ s(f) = (df /dt)~ 1 k on / g(f) . If one assumes that the time-dependence 
of the force is determined solely by the pulling, / = /i c ■ t, the loading rate (df jdt) = ji c 
reads /i c = k c ■ v. However, the situation becomes more involved if a soft linker such as 
a polymer, long nucleic acid, polysaccharide or polypeptide chain is coupled to the bond 
along the reaction coordinate[8j [TP] . In principle, the energy landscape can no longer be 
considered as one-dimensional, as the linker has to be treated explicitly. An example of a 
calculation for the escape from a single well has been presented in ref. [23] . In the present 
paper, we will discuss two approximative ways of treating the impact of a soft linker on 
the force spectrum. In the model of Evans and Ritchie [10] the change in the compliance 
of the system due to the presence of a soft linker gives rise to a modified loading rate: 

(df/dt)- 1 = ii- x [1 + C L (f)\ (9) 

with the bare cantilever loading rate fi c = k c -v and the force-dependent compliance of the 
linker C^. We will restrict ourselves to the worm-like chain (WLC) model for the description 
of a soft linker. In this case one has C WLC {f) = 2(3l p L c k c (l — z(f)) 3 /[l + 2(1 — z(f)) 3 ]. 
In this expression, z(f) is the force-dependent extension of the WLC-linker, i.e. z(f) is 
determined by the force law f(z) = (1/Af3l p )z/(1 — z) 2 (6 — 9z+4z 2 ). Furthermore, l p and L 
denote the persistence length and the contour length of the polymeric linker, respectively. 



c 
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Another way to treat the effect of a soft linker is the introduction of a potential of 
mean force where the linker coordinate is integrated out [8] . In this approach one considers 
the potential consisting of the sum of the potential given in eq.Q, Vo(q) + V pu ii(q,t), and 
the linker potential Vw^q — qi) with the linker coordinate q\. The potential of mean 
force is defined by the relation e-WMit) = Jdqie- m ^ q)+v ^ ll{q ' t)+Vlin ^ q ~ qi) \ In a harmonic 
approximation, this procedure results in a potential of the same form as given in eq.Q 
but with an effective force constant for the cantilever-linker system, 

k eS = (fc-* + kl 1 )- 1 (10) 

Here, k^ denotes the linker force constant. Therefore, in this approach one has to replace 
the bare cantilever force constant k c by a smaller k e s . We mention that for the WLC model 
the force constant of the linker is given by kwic = 3/(2/?/ p L c ). 

Note that in eq.(|8]) we already used the fact that the sum of the populations in well 
A and well B equals unity. The solution of the master equation is standard and can be 
performed in a semi-analytical manner. Knowledge of n(f) allows the computation of all 
quantities of interest. (If one is interested in the pull mode, one uses the natural initial 
condition n(f ) = 1 and in case of the relax mode, starting from a (high) force / max , this 
is replaced by n(/ max ) = 0.) 

In the standard description of dynamic force spectroscopy, one considers the rupture 
force distribution which is given by p(f) = —n'(f) with the prime denoting the derivative 
with respect to the force /. Note that this definition usually is employed in the absence of 
rebinding. However, we will use the same definition in case of finite k on , cf. ref.[3]. Hence, 
the rupture force distribution is given by: 

P(f) = ~i<S) = [««*(/) + V(/)] n(f) - K on (f) (11) 
Using this distribution one can calculate the moments 
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</ f n > = dff n p(f) = n dfr- l n{f) (12) 



o 



The most important ones, the mean rupture force (/) and the second moment of = ((/ r — 



(fr)) 2 ) will be discussed later. We will show that the mean rupture force defined by eq.(12), 
(f r ), behaves very similar to the maximum of p(f), f r ,m- 

As mentioned already in the Introduction, Seifert [T5] and also Li and Leckband[T5] 
considered a so-called dynamic strength, given by: 

F(f) = k c ([vt-q}) A (13) 

where ([v • t — q\)A = Z^ 1 J^dq [v • t — q) e~^ v ^ q ^ ■ n(f) and analogously for well B. Here, 
we will consider the averaged FE-curve 

F(f) = k c [([vt-q]) A + ([vt-q]) B ] (14) 
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because this can be obtained directly from simulations (and actual measurements near 
equilibrium [IT]). In a Gaussian approximation for the averages, cf. eq.(|B.l ), one finds: 



F(f) ~ [f - k c q A (f)Hf) + [f~ - "(/)) (15) 

where qx (/) = ikxq°x+f) / (kx+k c ) is the equilibrium position of well X under the influence 
of the external force. We will also briefly discuss the substantial differences between the 
local maximum of the averaged FE-curve, F m , and the maximum of the rupture force 
distribution, f r ^ m . 



III. Results and discussion 

A. Rupture force distributions and dynamic strength 

We start the discussion with the treatment of the Brownian motion in a potential of the 
form given in eq.Q. In Fig. 2a we show typical examples of FE-curves as obtained from 
Brownian dynamics simulations which were carried out employing a stochastic Runge- 
Kutta algorithm |23j. All parameters used in the simulation are given in the caption to 
Fig. la. These curves can be viewed as representative for experimentally observed FE-curves 
in the case of a hard linker. We have chosen a set of parameters for the potential closely 
related to what is known for the strength of a single hydrogen bond (Fig. la). Reversibility 
on experimentally relevant time scales was realized by placing well B energetically close to 
the transition state rendering the system irreversible only at high loading rates. 

It is evident that the rupture force and the rejoining force are stochastic quantities. Of 
course, whenever a rupture or rejoining event takes place, the system crosses the position 
of the transition state. This position is independent of force in our example, qr(f) = <?T' 
because we consider a steep barrier. This does not hold if a soft barrier is considered 
instead. We will later briefly discuss the case of a finite curvature of the potential at 
All further information about the system, such as rupture force distributions or averaged 
FE-curves, is based on an analysis of simulated FE-curves. 

In order to check the correctness of our definition of the rupture force distribution, 



eq.(ll), we compare the calculated distributions p(f) with the one obtained directly from 
the simulations (Fig.2b). It is evident that both methods yield the same results for high- 
loading rates. At smaller loading rates the situation becomes more complicated due to 
fluctuations in the force FE-curves rendering the correct determination of the rupture event 
difficult. A better way to validate the theory is hence to compare averaged FE-curves as 
demonstrated in Fig.2c. Perfect agreement between simulations and theory was found, as 
previously reported by Li and Leckband albeit we always use the factorization approxima- 
tion for the MFPTs, eq.Q. It also can be seen immediately that the mean rupture force 
increases as a function of loading rate in the pull mode and decreases in the relax mode, 
as has been found earlier [2]. Averaged FE-curves have been obtained by averaging 1000 
different individual FE-curves obtained by simulations, cf. Fig.2a. The dashed lines in 



Fig.2c are the Gaussian approximation to the average FE-curves, eq.(15). We found that 
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the agreement between the Gaussian approximation and the simulated curves decreases 
somewhat with increasing loading rate. However, for a loading rate as high as \i c = 3 • 10 5 
pN/s the maximum discrepancy is on the order of 5%. Given this satisfying agreement 
between simulation results and semi-analytical calculations, in the following we will solely 
discuss the results of the solution of the master equation ([8]), i.e. the distributions p(f). 

In Fig. 3a, we show rupture force distributions for various loading rates, both in the pull 
mode (upper panel) and the relax mode (lower panel). It is evident that the skewness of the 
distributions in the relax mode is just opposite to the one in the pull mode. Furthermore, 
we included the distributions for vanishing rebinding (dash-dotted lines). For loading rates 
\i c > 10 4 pN/s rebinding effects can safely be neglected for the parameters chosen. For 
small loading rates, however, rebinding is essential as it prevents the maximum of p(f) to 
reach zero instead of its equilibrium value. 

The rupture forces and rejoining forces are plotted versus loading rate in the upper 
panel of Fig. 3b. It is evident how the values from the pull and the relax mode converge to 
the equilibrium value for small loading rates. It is also obvious that the overall behavior 
of (/ r ) and f r ^ m is roughly the same. This does not hold for the mean rupture forces 
in the absence of rebinding, (f r )o (dotted lines). For small loading rates these behave 
qualitatively different, as equilibrium is absent. Only for large loading rates the effects of 
rebinding becomes negligible and (f r ) ~ (/ r )o- We mention that we calculated (f r )o in the 
pull mode by simply setting k on = in eq.pl). In Appendix C, it is briefly explained how 
the present off-rate is related to the one used in the calculation performed by Dudko et 
al.0. 

In the lower panel of Fig. 3b, we show the width of the rupture force distribution, 
a r = {(fr) ~ (fr) 2 ) 1 ^ 2 ■ The fact that a r is nearly constant for small [L c clearly reflects the 
influence of the finite rebinding probability. Only in the regime where rebinding becomes 
less important the width increases and for large \l c coincides with the corresponding one 
for vanishing rebinding (dotted lines). Therefore, in the interesting /x c -regime, the width 
looks similar to what is predicted by the Bell-model without rebinding[9]. Consequently, 
o r data should be treated with care. (The fact that the apparent slope changes for small /i c 
in the pull mode and for large [i c in relax mode is due to the fact that we only considered 
positive forces in our calculations similar to typical experimental situations.) 

As noted above, the rupture force distribution p(f) and the dynamic strength F(f) 
give rise to different maximum values for the 'rupture' force or the 'characteristic' force. 
This finding is a manifestation of the observation already made by Seifert [T8] that the 
two forces f r ^ m and F m are distinct quantities. We note that the differences between the 
averaged FE-curves and the dynamic strength depend on the loading rate. In Fig. 3c we 
plot the Gaussian approximation for F(f) versus the extension q for several loading rate 
and in the inset this is compared to F(f). In order to quantify the difference, we calculate 
various characteristic force values which are plotted in Fig.3d versus loading rate /i c . As 
indicated in Fig. 3c, in the pull mode the F m -values are smaller than F m and these in turn 
are smaller than f r>m . When the data are shifted to start at the same low-/i c limit, one 
can see the different apparent slopes of the various characteristic force values. For large 
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[Ac, the difference in slope between / nm and F m is on the order of 25%. In the relax mode, 
the differences between the various values are much less outspoken, mainly due to their 
smaller overall variation. From an experimental point of view, systems close to equilibrium 
displaying FE-curves fluctuating between state A and B are better analyzed in terms of the 
averaged FE-curves instead of the rupture force distributions p(f) by recording merely the 
last rupture or rejoining event, as already mentioned by Evans and RitchiepQ. Of course, 
the differences between the characteristic forces F m and / nm have to be taken into account. 

We close this section with noting that the differences in the 'rupture' and 'rejoining' forces 
determined from averaged FE-curves and from rupture and rejoining force distributions 
have their origin in the intrinsic nonlinearity of all force-dependent quantities. This also 
holds in equilibrium, see Fig. 3d. We mention that it is not simply related to the non- 
gaussian shape of p(f), as can be shown by calculating F(f) with an assumed Gaussian 
n(f) (not shown). 

In the following we discuss the dependence of experimentally relevant properties such as 
rupture and rejoining forces on kinetic parameters. We mainly focus on such quantities that 
have not already been discussed earlier by Li and Leckband for the dynamic strength |19j. 
However, in the following, we will use Kramers rates for a general bistable potential in order 
to be able to also discuss the influence of a soft barrier. We will focus on the distribution 
of rupture forces and quantities derived from p(f). 



B: Variation of kinetic parameters 

In order to be able to include a discussion of the effect of a soft barrier and not to be 
restricted to the quite special form of the potential discussed above, in this section we 
consider a more general form of a double-well potential and apply Kramers theory to the 
calculation of the transition rates. Without going into details, we give the results for the 



off- and on-rates in Appendix A, eq.(A.5). The parameters can be chosen in such a way 
that both, the shape of the rupture force distributions and the mean values coincide with 
those obtained for the cusp-like potential discussed above. This is demonstrated in Fig.4a, 
where (/ r ) is plotted versus /i c . Also included in that figure is the result of a calculation 



with rates according to the Bell model, cf. eq.(A.8)[6J. The discrepancy between the latter 
and the local harmonic or the cusp model is similar to what is observed in case of a pure 
escape model [9|. 

Next, we consider the parameters relevant for the characterization of the transition 
state in the local harmonic approximation, namely the position and the curvature kr- 
The mean rupture forces (f r ) are plotted versus \i c in Fig. 4b for different values of these 
parameters. In the upper panel q% is varied and it is obvious that the effect is just what one 
expects from the Bell model, namely the (/ r )(log(// c ))-curves become steeper for smaller 
q\. We note that the discrepancy between the local harmonic approximation and the Bell 
model strongly depends on the distance of the well to the transition state. In order to 



show this, we included results from calculations using the Bell expressions, eq.(A.8) as the 
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dotted lines in Fig. 4b. It is evident that the discrepancy diminishes with increasing q^. The 
reason for this behavior is easily understood from the derivation of the Bell expressions for 



the transition rates as the limit of large critical force f c of the Kramers rate, cf. eq.(A.8), 
because f c is directly proportional to q°r. Explicitly, one has for large fi c , where rebinding 
is negligible, cf. e.g. [9]: 

(/r)pull - and {frlrelax - ~ (16) 

QTA WBT 

with Zp = e'^QrA/k^g, z r = e~ 7 QBT/k® Q , 7 = 0.577215... being the Euler constant. Of 
course, in the relax mode the effect is very small for the parameters chosen. 

This situation changes when a variation in curvature is considered (Fig. 4b, lower panel). 
In the Bell model, the on- and off-rates and therefore also the rupture and rejoining forces 
are independent of the value of kx- Here, we find an increasing apparent slope at large 
loading rates for decreasing kx because we are using the on- and off-rates in the Kramers 



approximation, eq.(A.5). At small loading rates, however, (f r ) is hardly affected at all. 
In the relax mode, no changes in the rejoining force are found for the present parameters. 
This behavior can be understood qualitatively from the expressions for the scales of the 



relevant forces given in eq.(A.6). For the off-rate this scale is given by f c ~ k!+k T $ anc ^ 
thus directly proportional to far- In the relax mode, the relevant rate is the on-rate and 
here one has for the scale f r ~ (k B + k c )(q B — q T ) independent of k?- (Note that this 
argument holds for soft B- wells only.) 

The dependences of the characteristic scales f c and f r on the various force constants 
also allows to qualitatively understand the dependence of the rupture and rejoining forces 
on the cantilever spring constant k c for large \i c . Again assuming ik^ ~ k T ) 3> {k B ~ k c ) 
we have almost no dependence of k oS (f) on k c whereas k on (f) changes strongly with k c . 
Thus, for large /i c we expect a very weak k c dependence of (f r ) in the pull mode and a much 
stronger one in the relax mode. This is just what is observed in the calculation as shown 
in Fig.4c. In the upper panel one can see that for large /i c the behavior just discussed is 
found, whereas for small fi c there is strong dependence on k c , also in the pull mode. In 
the lower panel, the dependence of the mean rupture force in equilibrium (/i c = 0) on k c 
is displayed. For comparison, also the results for the Bell model are shown as the dotted 



lines. In case of this model the rates are given by eq.(A.8) and the mean rupture force can 



be calculated analytically (cf. eq.(12)): 



(f r r= 4fn«K/)=/ 7 =—111(1 + ^(0)) (17) 

Jo Jo l + K eq L (f) Q BA 

where K eq (f) = k on (f)/k oS (f). Similarly, one finds f^ m = In {K cq (0))/Q B A- It is important 
to note that K eq (Q) for the parameters chosen depends on k c . The same holds for Qba = 
(iks/iks + k c )q%. The upper dotted line has been calculated using this expression and the 
lower one with the approximation Qba = flq B - (If one uses t ne k c = 0- value for K eq (0) 
along with Qba/ (3 = 1 nm, one finds (f r ) eq — 26.3 pN for the parameters chosen.) One can 
readily see that the distance between the two minima can be extracted from (/ r ) eq data if 
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these are extrapolated to k c — > 0. In this case all three curves yield a distance between q A 
and q Q B of about 1 nm within 20% accuracy. Thus, if one has an estimate of the equilibrium 
constant K eq (0) from other kinetic measurements, it is possible to get information about 
the position of the minima q° A and q% from the equilibrium mean rupture forces. On the 
other hand, as the mean rupture forces in the large /z c -limit allow estimates of Qta and 



Qbt, cf. eq.(16), one can use the fc c -dependence of the equilibrium mean rupture force in 



order to determine K eq . 



C: Impact of a soft linker 



As has been mentioned in the preceeding section, one frequently is confronted with a 
situation where the influence of a soft linker cannot be neglected. In the model considered 
by Evans and Ritchie [TU], for vanishing rebinding the main effect of a soft linker is to 
broaden the rupture force distributions somewhat and to decrease the mean rupture forces 
due to the smaller effective force constant of the cantilever-linker system. We use eq.([9]) 
in order to treat a soft WLC-linker in conjunction with the master equation, eq.(|8]), for 
finite rebinding rates. In Fig. 5a, we show the influence of polymeric linkers of varying 
contour length L c on the rupture force distributions. It is seen that the effect of the 
smaller compliance of the linker is very similar to what has been found in the absence of 
rebinding in the pull mode, i.e. a shift of the peak maximum towards smaller forces. In 
the relax mode just the opposite shift is found, as expected. Therefore, the mean rupture 
forces determined in the two modes shift together and the equilibrium force is reached at 
a somewhat higher loading rate than in case of a hard linker. Note the distortion of the 
p(f) curves indicating more weight at lower (pull) or higher (relax) forces. It thus appears 
that a soft linker has some impact on the skewness of the rupture force distributions. 

The shift in the mean rupture forces is shown in Fig. 5b, where we plot (/ r ) versus 
\i c for a hard linker (L c = 0) and a WLC-linker with a contour length of 50 nm and a 
persistence length of 0.4 nm. Also the fact that the rupture force in equilibrium, (/ r ) eq , 
is independent of the linker is obvious. This fact can easily be understood from the equi- 



librium solution of the master equation, cf. eq.(17), n(f) ccl (f) = [1 + K eq 1 (f)] 1 with 
K cq (f) = K on (f)/K oS (f) = k on (f)/k oS (f) independent of (df/dt) and thus of C L (f). There- 
fore, the rupture force (/ r ) eq is unaffected by the influence of a soft linker as observed for 
the B-S transition of DNA[25]. 

However, we note that this is not a generally valid result. If we treat the influence of 
a soft linker in the manner discussed by Hummer and Szabo[8], we have to use a small 



effective force constant given by eq.(10) instead of the bare cantilever force constant. The 
dotted lines in Fig. 5b are obtained this way and one clearly sees a behavior similar to that 
discussed in the context of the fc c -dependence of (f r ) eq , cf. Fig.4c. 

From this discussion we conclude that measurements of the equilibrium rupture force 
can in favourable cases be used in order to improve the understanding of the impact of soft 
linkers on the force spectrum. 
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D: Bond heterogeneities 

It has been shown that in some cases the experimentally observed skewness of the rupture 
force distributions is opposite to what is expected on the basis of model calculations, 
independent of the model used [131 EH]. The explanation put forward by Raible et al. is 
that the parameters entering the off-rate are distributed. In particular, if the off-rate 
is written as k Q s(f) = k® s e QrA f , a distribution of the distance to the transition state, 
Qta-i is considered and this is compatible with the experimental findings for a number of 
examples [26J. 

If one considers a finite rebinding probability with a rate k on (f) = k® n e~® BT f one has 
two parameters Qta and Qbt that can be considered as distributed. We have performed 
calculations of the rupture force distributions and the mean rupture forces with distribu- 
tions of both parameters. Without showing results here, we mention that the effect is 
very similar to what has been observed in the case without rebinding. The widths of the 
rupture force distributions are more strongly affected by bond heterogeneities in the pull 
mode than in the relax mode. The mean rupture forces are hardly changed at all. 

III. Conclusions 

The model calculations presented in the present paper demonstrate the importance of finite 
rebinding probabilities in the interpretation of dynamic force spectroscopy experiments. We 
have shown that the usual definition of the rupture force distribution also is meaningful if 
rebinding cannot be neglected. Therefore, one can discuss the so-called force spectrum in 
the same way as it is usually done when one is concerned with irreversible bond breaking 
events. There is, however, a difference to the dynamic strength, a bulk property, that 
has been considered in earlier work on reversible kinetics in the context of dynamic force 
spectroscopy. From an experimental point of view, however, near equilibrium fluctuations 
in the individual FE-curves can prevent correct data aquisition for rupture and rejoining 
force histograms. In such cases it may be more practicable to record averaged FE-curves 
as shown by Bornschlogl and Rief|17j. 

One effect of finite on-rates regards the interpretation of the width of the rupture force 
distributions. It has been shown earlier [9] that this width increases with increasing loading 
rate without rebinding. Here, we find that due to finite on-rates the width is almost 
constant over a larger range of loading rates and only for large fi c starts to increase in the 
same way as without rebinding. Therefore, in order to discriminate among various models 
for the transition rates one has to carefully cover a large range of loading rates. 

We expect that it is not easy to discriminate between the Bell model and more complex 
scenarios for distances to the transition state typical for proteins. This is due to the fact 
that the Bell-like rates are obtained from those for a local harmonic approximation in 
the limit of large critical forces. These critical forces themselves are proportional to the 
distance to the transition state. For larger distance the force spectrum shows less curvature, 
cf. Fig. 4b. 
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It has been shown that the determination of the equilibrium rupture force allows the 
extraction of the equilibrium constant K eq at zero load. This quantitiy, however, depends 
on the force constant of the cantilever and should therefore be extrapolated to k c — > in 
order to be able to compare it directly to other determinations like bulk measurements. 
Small cantilever force constants are generally recommended alos to minimize the impact 
of polymeric linkers and their variation in contour length on the outcome of a force exper- 
iment. 

We investigated the impact of a soft linker on the rupture force distributions and 
the mean rupture force. We found that in principle it should be possible from detailed 
determinations of the equilibrium rupture force to gain some deeper understanding of the 
compliant behavior of the composite cantilever-linker system. 

In conclusion, we have shown that dynamic force spectroscopy can yield insightful 
results on the nature of system showing reversible bond-breaking and that the equilibrium 
regime of small loading rates is particularly interesting in this respect. 



Appendix A: Transition rates 

MFPTs for the harmonic cusp-like potential 

In this Appendix, the explicit expressions for the transition rates for the cusp-like potential 
and the Kramers rates in a general potential are given. For a potential built from two 
harmonic wells with a sharp intersection on finds from eq.([5]): 

r%T dUC \f) = D- l Z A I AT and r^ c \f) = D- l Z B I BT 



with 



Zx(f) 



-PVxU) 



• 27T \ 1/2 



71 



1/2 



2(3k 



xc , 



erfc{A XT (/)} 



X = A, B 



Here, we used the abbreviations kxc = kx + k c 



VxU) = Vx-\kxcqxU) 2 + \kx(qx? with qx(f) = k? c (kx& + f) 



and 



A 



'Pk 



AC 



1/2 



AT 



fer(/) ~ q A (f)} 



A 



'(3k 



BC 



1/2 



(A.l) 



(A.2) 



BT 



[qnif) ~ q T (f)} 



Furthermore, we defined 



Ixr(f) 



( 2 \ 1 ^ 2 r x 

— e^ Vx{f) d(A XT (f)) with d(x)= dze z 

\pkxc J Jo 



2x 



(A.3) 



Using eq.(A.l) yields results for the populations of the potential wells as a function of 
force (time) that are indistinguishable from calculations with the exact expression for the 
MFPTs. 
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The Kramers ap proxi mation for the MFPTs is obtained from eq.( A.l ) if one keeps only 

I I _2 

the first term in eq.(A.2) and using d(x) 



2x 



Kramers 
'XT 



7T 



exp (A 



It) 



j3Dk 



xc , 



A 



XT 



(A.4) 



Transition rates in a local harmonic approximation 

Next, we give the expressions for the Kramers rates for off- and on-rates in a general double 
well potential [21, 22J. Approximating the minima as well as the maxium by a parabola 
with curvatures Jsa, k B for the minima and fcr for the maximum, one finds: 

KM = k^^-WM and U/) = Ce- WflT/[l+i///p] (A-5) 

Here, Qta = Qt - Qa, Qbt = Qb - Qt with Q x = kx/kxcQx for X = A, B and 
Qt = hr/ (kr — k c )q T . Furthermore, the scale of the forces are set by 

{k T -k c )k AC ^ , , {.k T -k c )k BC 

fc = — : —j Qta and f r = — — Qbt (A.6) 

kA + k T k B + k T 

and the rates in the absence of any force are given by: 



k° - n 



^{k T - k c )k A c 
2^ 



e -l3(V T -V A ) and k C 



D 



\J(k T - k c )k B c 
2^ ' 



-P(V T -V B ) 



(A.7) 



Finally, it easy to show from eq.(A.5) that in the limit of large f c and f r the rates of the 
Bell model 

S'(/) = Ce WH/ and *£>(/) = k on e'^f (A.8) 

are obtained. Here, Qta and Qbt in the limit k c — > are the respective distances to the 
transition state. 



Appendix B: Averaged force versus extension curves 

Here, we derive the expressions for the averaged FE-curves from the rupture force distri- 
bution. In a Gaussian approximation, the expectation values (q(f)) are given by 

W))x = Z-x l I dqe-M'fiq = q x (f) = **f±l (B.l) 
Jx k X c 

with X = A or X = B depending on the potential well considered. In this approximation 
the increase in the force as a function of the extention is simply given by [/ — k c qx(f)], 
which of course is the same as the force measured for instance by an AFM cantilever as 
long as the bond considered is intact. With this, the overall force is given by: 

- k c q A (f) for f < f r , . 

-k c q B (f) for/>/ r ^ 
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This means, one can write 



F(f, f r ) = [f~ k c q A (f)Wr ~f) + [f- kcMf)Mf - fr) (B.3) 

where Q(x) denotes the Haeviside step function, Q(x) — 1 if x > 1 and = if x < 0. 
Here, f r is the rupture force and eq.(B.3) is just an expression of the fact that after rupture 
the system is in the right hand minimum of the potential. The average over the rupture 
forces is given by the rupture force distribution, i.e. 

(F(f)) = / df r F(f, f r )p(fr) = [f- k c q A (f)Hf) + [f- k c q B (f)](l - n(f)) (B.4) 
Jo 

which follows immediately from the definition of Q(x) and the fact that p(f) = — -^n(f). 



Appendix C: Escape from a harmonic well 

Here, we discuss the transition from a double well potential to a model of escaping a single 
harmonic well as discussed by Dudko et al.[9j. In the product approximation of eq.([5]) one 
has for the off-rate, cf. eq.([6]): 

Ki(f) = r AT (f) + r BT (f)Z A (f)/Z B (f) - D- l Z A [I at + Ibt] (CI) 

If the potential is given by the sum of eq.Q and the pulling potential according to eq.Q, 
the limiting case of a single well adjacent to an unbound state is provided by the limit 
Vb — > (— oo). If one chooses a fixed value for the force constant ks, the condition of the 
crossing point of the two parabola yields q% = ([2(Vt — Vb) — ^(gy) 2 ]/^) 1 ' 2 where q? 
is defined by Vr = l/2k A c{qj) 2 ■ For large negative Vb one then has V B — —l/2k B (q B ) 2 . 
If one now uses eqns.( [Al2|A.3 ) one finds that 



I BT ^ {(3k B q B Y l and thus kJ{f)~r AT (f) (C.2) 
The resulting rate for escape from a harmonic well, k s (f), can be cast into the form[9j 



k a (f) = k a (0)e^ 1 -^ with fc,(0) = D i(3k fl q° T e-^ (C.3) 

V 27T 



where the abbreviation e = 1 — //(A;ac , ?t) nas Deen used and Vp is given above. 
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Figure captions 



Fig.l (color online) 

a: A sketch of a double well potential as used in the calculations, the explicit 
parameters are chosen as follows: The cantilever force constant was set to k c = 30 
pN/nm, k A = 1000.0 pN/nm and k B = 48.49 pN/nm; the diffusion constant is 
D = 2000.0 nm 2 /s; the potential heights are chosen as V A = 0, V B = 8.0/(3 = 33.12 
pNnm and from this one obtains the barrier height Vr = {k A + k c )(qj,) 2 /2 = 46.35 
pNnm; the force-free equilibrium positions are q A = 0, q?> = 0.3 nm and q% = 



q° T + sJ[2{V T - V B ) - kM) 2 ]/k B = 1.0 nm; 

b: Kinetic rates k on / ff(f) vs. the force / in various approximations. The inset shows 
the ratio of the exact off-rate to the product approximation ^ and the Kramers 
approximation, indicating the quality of the former. 

Fig. 2 (color online) 

a: Two examples of FE-curves as obtained from Brownian dynamics simulations. 
The potential is of the shape shown in Fig. la and the explicit parameters are the 
same as in Fig. la. 

b: Rupture force distributions for the pull mode, p(f), and the relax mode, p Te i(f), 
as obtained from Brownian dynamics simulations (histograms) and from the solution 
of the master equation, eq.([8]). 

c: Dynamic strength from Brownian dynamics simulations (full lines) and from the 



Gaussian approximation, eq.(15) (dashed lines). Upper panel: fj, c = 300 pN/s, lower 
panel: \i c = 3000 pN/s. 

Fig. 3 (color online) 

a: Rupture force distributions for the pull mode (p(f), upper panel) and the relax 
mode (prei(/), lower panel) versus force / = \i c -t for /i c — 10 2 , 10 3 , 10 4 , 10 5 pN/s (for 
p(f) from left to right and for p Te \(f) from right to left). The remaining parameters 
are the same as in Fig. la. Also shown are the distributions for the case of vanishing 
rebinding, i.e. fc° ff = for the pull mode and k® n = for the relax mode as the dash- 
dotted lines. The equilibrium rupture force distribution p cq (f) = p^iif) is shown as 
the dotted line. 

b: Rupture and rejoining forces (upper panel) and variances (lower panel) as a func- 



tion of the loading rate. The full lines are the mean rupture force according to eq. ( 12 ) , 
(f r ), and the dash-dotted lines represent the values obtained from the maximum of 
the rupture and rejoining force distribution, f rm . The dotted lines are the mean 
rupture forces without rebinding, (f r )o- In the lower panel a r = J (/ 2 ) — (f r ) 2 is 
shown including the case of vanishing rebinding (dotted lines). 



c: Averaged FE-curves calculated in the Gaussian approximation, eq.(15), for pull 
mode (full lines) and relax mode (dash-dotted lines), \i c = 10 2 , 10 3 , 10 4 pN/s. The 
maxima increase for the pull mode and decrease for relax mode. Also shown as the 
dotted line is the equilibrium curve, /i c = 0. The two straight lines are the factors 
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(f-k c qx(f)) 



X = A, B. The inset shows the comparison to the dynamic strength 



F(f), cf. eq.(13), for the pull mode (dashed lines). 

d: Rupture forces / rjTn (full lines), F m (dashed lines) and F m (dotted lines) as a 
function of the loading rate. The inset shows the data in the pull mode with F m and 
F m shifted in order to match f r ^ m at small \i c . 

Fig. 4 (color online) 

a: Rupture force (f r ) versus loading rate in the local harmonic approximation 



(Kramers rates, eq.(A.5), full lines), the Bell model (eq.(A.8), dashed lines) and for 
the cusp potential discussed in subsection A. The parameters for the Kramers rates in 
the local harmonic approximation are: the local curvatures are kA = 1300.0 pN/nm, 
ks — 48.49 pN/nm and kx = 10 4 pN/nm; the diffusion constant is D = 1442.75 
nm 2 s; the potential heights are Va — 0, Vb — 33.12 pNnm and Vp = 46.35 pNnm; 
the force-free equilibrium positions are q A = 0, q T = 0.3 nm and q B = 1.0 nm; note 
that Va, Vb, Vt, q% 9t = an d q% are the same as in Fig. la; the cantilever force 
constant was set to k c = 30 pN/nm. 

b: Effect of variation of the parameters characterizing the transition state on (f r ): 
upper panel: Variation of the position qj. for force constant kx = 10 3 pN/nm. Also 
included as the dotted lines are the results for Bell rates in the pull mode for large 
\i c for q T = 0.3, 0.4 and 0.5 nm. 

lower panel: Variation of the curvature kx for constant position q T = 0.4 nm. The 
remaining parameters are the same as in Fig.4a. 

c: Mean rupture force for varying cantilever force constant k c . Upper panel: (f r ) ver- 
sus \i c for different k c . Lower panel: (f r ) eq versus k c . For comparison, also the results 



for the Bell model, eq.(17), are shown as the dotted lines: Qba = Pk B /(kB + k c )q B 



(upper line), Qba — Pq% (lower line). All other parameters are the same as in Fig.4a. 
Fig. 5 (color online) 

a: Rupture force distributions in the pull mode and the relax mode in the presence 
of a soft WLC linker for different loading rates fi c . The persistence length is fixed to 
be l p = 0.4 nm and the cantilever force constant used is k c = 5 pN/nm. The other 
parameters are the same as in Fig.4a. The different lines correspond to different 
contour lengths L c as indicated in the uppermost panel. Here, L c = means that no 
soft linker is included, cf. eq.(|9]) and the expression for Cwlc- 

b: (fr) eq versus fi c for contour lengths L c = (full lines) and L c = 50 nm (dashed 
lines). The other parameters are the same as in (a). Also included as the dotted 
line is the result from a calculation with an effective cantilever force constant k e g = 
(/c" 1 + fc^/LcO^E! where k WLC = 3/(2/3l p L c ) ~ 0.3105 pN/nm is the force constant 
of the WLC-linker. 
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